To APPEAR IN ApJ Letters 

Preprint typeset using I4TgX style emulateapj v. 1 0/09/06 



CORRECTING STELLAR OSCILLATION FREQUENCIES FOR NEAR-SURFACE EFFECTS 



00 

o 
o 

(N 



To appear in ApJ Letters 

ABSTRACT 

In helioseismology, there is a well-known offset between observed and computed oscillation frequencies. 
This offset is known to arise from improper modeling of the near-surface layers of the Sun, and a similar effect 
must occur for models of other stars. Such an effect impedes progress in asteroseismology, which involves 
comparing observed oscillation frequencies with those calculated from theoretical models. Here, we use data 
for the Sun to derive an empirical correction for the near-surface offset, which we then apply three other stars 
{a Cen A, a Cen B and (3 Hyi). The method appears to give good results, in particular providing an accurate 
estimate of the mean density of each star. 

Subject headings: stars: individual (/3 Hyi, a Cen A, a Cen B) — stars: oscillations — Sun: heUoseismology 
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1. INTRODUCTION 

Both helio- and asteroseismology involve comparing ob- 
served oscillation frequencies with those calculated from the- 
oretical models. However, for the Sun there is a long- 
standing systematic offset between observed and computed 
frequencies that is know n to arise from improper modeling 
of the near-surface layers (Christensen-Dalsgaard et alj 
Dziembowski et al. 1988; Christensen-Dalsaaard et al 



scale (see also lChristensen-Dalsgaard & Goughlll980l) : 



1988 



1996 



Christensen -Dalsgaard & Thompson 1997). This offset is in- 
dependent of the angular degree of the mode (/) and increases 
with frequency. A similar offset must occur for models of 
other stars, and should be taken into account whenever ob- 
servations of stellar oscillations are compared with theory. In 
this Letter we use data for the Sun to derive an empirical cor- 
rection for these near-surface effects and show how to apply 
this to other stars. The method appears to give good results, in 
particular providing an accurate estimate of the mean density 
of each star. 



2. METHOD 

The p-mode oscillations in solar-type stars for a given an- 
gular degree I are approximately equally spaced in frequency, 
with a separation of Af (the so-called large separation; see 
IChristensen-Dalsgaard 2004 for a review of the theory of 
solar-like oscillations). Since the offset from incorrect mod- 
elling of the near-surface layers is independent of /, we can 
derive the correction by considering only the radial modes 
(/ = 0) and then apply it to all modes. 

Suppose we have a set of observed frequencies for radial 
modes, i'obs(n), where n is the radial order Suppose also that 
Z4,est(n) are the frequencies from the best model, by which we 
mean the one that best describes the parameters and internal 
structure of the star, but which still fails to model correctly the 
surface layers. 

For the Sun, the difference between observed and best 
model frequencies turns out to be well fitted by a power law, 
which has the convenient property of being free of a frequency 
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fobs(n)-t'best(«) = fl 



'^obs(«) 



(1) 



where t'o is a suitably chosen reference frequency, and a and 
b are parameters to be determined. 

Suppose further that we have calculated a reference model 
that has frequencies i'ref(n) and is close to the best model. 
From homology scaling it then follows that, to a good ap- 
proximation, 

TAxJn) = ny,st(n), (2) 

where the scaling factor r is related to the mean densities pbest 
and pref of the best and reference models by 

Pbest = r Pref- (3) 

Given a determination of r, equation (O provides our best es- 
timate of the mean density of the star 

Substituting equation ^ into equation ([T]) gives 



i^ohs(n)-ny,sf(n) = a 



1^0 



(4) 



and differentiating with respect to n gives 



A ^ \ A r \ A/^^obs(n)\ Avobs(n) 
Ai^ob,s(n) - '"Ai'ref(n) = ab . 

Combining and rearranging these last two equations gives 

fref(n) Al/,.ef(n) 



r = {b-l)[b 



l/obs(«) Aj/obs(«) 



(5) 



(6) 



and 



Al'ref(n) 



t'ref(n) 



(7) 



Ai'obs(n) V V ^obs(n) 

If we know b then we can calculate r using equation (|6]l, or 
vice versa using equation (Q. We can then obtain a using 
equation ^. 

We now show how to apply the above method to a set of 
observed and calculated frequencies. Suppose we have fre- 
quencies for A^ radial modes with orders ni,n2,... ,11^ (not 
necessarily consecutive). We use these to calculate the four 
terms needed to evaluate equation (|6]l or (|7]l- For i'obs(«) and 
i/ief(n), we simply use the means of the given sets of frequen- 
cies, which we denote by (i'obs(n)} and (i'ref(n))- 
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To estimate the large separations, we calculate the slope of 
a linear least-squares fit to the given frequencies (as a function 
of n): 

J2ti ('^obsC";)- (i'obs(n))) («;- («)) 



{Ai^obs(n)> = ■ 



(Az/,-ef(«)) = 



J2tl ('^i-eK";)- (t^ref(»))) (»/-(»)) 

Em(«|-(«)) 



(8) 



(9) 



We must then assume a value for either b or r, and use equa- 
tion (|6]l or (|7]) to estimate the other. Finally, the value of a is 
found from equation (|4]l, as follows: 



a = 






f 



(10) 



We now proceed to apply this method to the Sun in order to 
measure b {^, and then adopt this value of b for other stars 
(®. 

3. APPLICATION TO THE SUN 

For the Sun we too k Model S of 
IChristensen-Dalsgaard et all (Il996l) . as listed in the first 
row of Table[Tl We assumed this to be the "best" solar model, 
in the sense defined above, which means we can set r = 1 
(see equation |2]i. For the observed solar frequencies, we 
used those measured by iLazrek et alJ (Il997h with the GOLF 
instrument on the SOHO spacecraft. 

We followed the procedure described above, choosing t^o = 
3100/iHz and setting r = 1, and using the data to measure b 
and a. We have chosen to use the nine modes centred at the 
peak of the oscillation power, from which we obtained the re- 
sults shown in the first line of Table|2] and a value of b = 4.90. 
The differences between observed and Model S frequencies 
are plotted as the squares in Fig.[Tl and the solid curve is the 
function given by equation ([T]i- 

The above fit was made for the strongest 9 radial modes 
in the Sun. We repeated the analysis for different numbers of 
modes (all values from 7 to 13), and found the derived value of 
b to range from 4.4 to 5.25. Clearly, the frequency differences 
do not exactly follow a power law, and so the exponent in the 
power-law fit depends substantially on the frequency range. 
Importantly, the value for a varies by less than 0. 1 /xHz in all 
cases. 

In addition to Model S, we also considered models de- 
noted S~ and S""" from the same evolution sequence as 
Model S, but with ages of 2.25 and 7.44 Gyr, respectively. 
The parameters of these models are given in Table[T] We kept 
b fixed at the value found for the "best" model (4.90) and used 
equation (|6]l to estimate r. The results are shown in Table |2] 
and also in Fig. [T] Importantly, the derived density of the 
Sun (r^pref) is correct for both these calculations, despite the 
very different densities of the models themselves, giving us 
confidence that the method has been successful. Even more 
importantly, the derived density of the Sun is completely in- 
sensitive to the choice ofb. This reflects the fact that the value 
of r obtained by fitting to equation|4]is not sensitive to the ex- 
act form of the function on the right-hand side, so long as that 
function tends to zero with decreasing frequency. 

4. APPLICATION TO a Cen A, a Cen B AND (3 Hyi 

We have considered three stars for which observations and 
models have been published. For a Cen A we took ob- 
served frequencies (radial modes only) from four sources: 



lBouchy_& Carried (|2002l), iBe dding et all (|2004|) . iBazot et al.l 
(2007) and Fletcher et al. (20061). This gave a set of 33 mea- 
sured frequencies for 1 1 orders, all of which were given equal 
weight in the fitting proc ess. For a Cen B th e observed fre- 
quencies were t aken fromlKieldsen et al.l (2005) and those for 
/3 Hyi from Bedding et al.l (l2007h . In each case, we used all 
the detected 1 = frequencies. 

The models that we have used for these stars are listed in 
Ta ble [H These includ e p ubhshed models of a C en A and B 
by Morel et al. (2000) and lThevenin et al." (2002), of a CenB 
by Tang et al. (2008), and of f3 Hyi by Fe rnandes & Monteirol 
(2003). In addition to published models, we have consid- 
ered seve ral computed with the Aarhus stellar evolution code 
(ASTEC. IChristensen-Dalsgaardll2008h . Models A and B ai-e 
models of the a Cen system, computed with essentially the 
same physics as Model S and fitted (T. T eixeira et al. 200 8, in 
preparation) to the ob servations of Bedd ing et al.l (|2004|) and 
iKieldsen et aTl (l2005l). For /3 Hyi , Model H matches the pa- 
rameters reported by iNorth et aJ] (|2007) and was computed 
with similar physics, but neglecting diffusion and settling. 
Finally, Models H" and H"*" are from the same evolution se- 
quence as Model H, but with substantially different ages and 
densities, which bracket those of Model H. 

For each model we used the value of b found for the Sun, 
and used equation ^ to estimate r. The results are shown 
in Table |2l and also in Figs. |2]-H] Note that there is consid- 
erable scatter in the observed frequencies of a Cen A and B 
due to the relatively short span of the observations relative to 
the mode lifetime. Taking this into account, it is once again 
encouraging to see that the power law (with a single value of 
b) provides a good fit to the frequency differences, and that 
for each star there is good agreement between the densities 
derived from the different models. 

To estimate uncertainties, we have repeated the fits for the 
same range of b values considered in ^(4.4-5.25) and again 
found that there is no effect on the calculated density, to the 
precision quoted in Table|2] Over this range, the change in a is 
less than 1 fiHz for /3 Hyi, less than 0.4 /iHz for a Cen A and 
less than 0.7 /xHz for a Cen B. The changes in the correction 
terms plotted in Figs.[2}Slare comparable or slightly bigger 

For each star, we can identify the model that is closest to be- 
ing the "best" model as the one having r closest to unity (note 
this is not necessarily the model with the smallest near-surface 
offset). These give our best estimate of the stellar density for 
each s tar. 

Nor th et al] (l2007h recently used interferometry to measure 
the angular diameter of f3 Hyi to be 2.257 ± 0.019 mas. They 
combined this with the parallax (from Hipparcos) and the 
mean density (from asteroseismology) to determine the ra- 
dius and mass of the star Using our best estimate of the 
mean density of /? Hyi (0.258 ±0.001 gcm"^ ') and the re- 
vised iZ/pparcoj' parallax (134.070 ±0.110 mas; lvanLeeuwenl 
l2007h . we derive slightly updated values, finding a radius of 
1.809 ±0.015/^0 (0.85%) and a mass of 1.085 ±0.028Mo 
(2.6%). 

5. DISCUSSION AND CONCLUSIONS 

The method outlined here for correcting near-surface ef- 
fects can be applied to model frequencies before they are com- 
pared with observations. As in the case of the Sun, we expect 
that the correction is independent of degree at a given fre- 
quency, for low-degree acoustic modes, and thus the correc- 
tion determined from radial modes can be applied to all such 
modes. 
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There is, however, an important exception. In evolved 
stars, mixed modes may be observed that have the character 
of gravity modes in the deep interior of the star. Observa- 
tional evidence for mixed modes has been foun d in ?/ Boo 
pCieldsen et al.' 19951 EQQI [Carrier et alJ 12005) and (3 Hyi 
(iBedding et al. 2007) . Owing to the larger amplitude of these 
modes in the stellar interior and hence their higher inertia, 
their frequencies are less affected by the surface effects. Tech- 
niques need to be developed to take this into account in the ap- 
plication of the surface correction before the frequencies are 
analyzed. 

The use of a single power law is made plausible by the 
knowledge that the offsets we are modelling are caused 
by the properties of the near-surface layers and hence pre- 
sumably depend only on surface gravity, effective temper- 
ature and composition, and not on the details of the in- 
ternal properties of the star This may be particularly 
true of the exponent which, in the simple analysis by 
IChristensen-Dalsgaard & Goughl ( 119801) . is determined by an 
effective average polytropic index (i.e., the relation between 
pressure and density) in the near-surface layers. 

It is important to note that, owing to its strong frequency 
dependence, the offset also affects the large frequency sepa- 
ration Al/, as is indeed implicit in equation (|5]l. For example. 



although Model S is one of the best available models of the 
Sun, it has a large separation that is 1 /iHz greater than the ob- 
served value. Thus, attempting to fit to stellar models based on 
Aiy will introduce systematic errors, unless the corresponding 
correction is applied. 

We have shown how to identify the model for each star that 
is closest to the "best" model, by requiring that r be as close 
as possible to unity. This gives us an extremely accurate es- 
timate of the mean stellar density. However, it is important 
to point out that, while a model with r close to unity gives a 
good match to the overall structure of the star, it does not nec- 
essarily reproduce the structure of the core or give a reliable 
estimate of the stellar age. Determination of those properties 
requires taking into account the frequencies of the non-radial 
modes (including the small frequency separations). That is 
the next stage of model fitting, which can be done after the 
surface correction has been applied to all modes. 
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TABLE 1 
Details of theoretical models 



Model ID 


Star 




Mass 


Radius 


Luminosity 


Age (Gyr) Reference 


Model S 
Model S- 


Sun 


1.000 
1.000 


1.000 
0.930 


1.000 
1.286 


4.52 
2.25 


Christensen-Dalssaard et al. (1996) 
this paper 


Model S+ 




1.000 


1.128 


0.832 


7.44 


this paper 


FM2003 
Model H 


/3Hyi 
/3Hyi 


1.100 
1.080 


1.899 
1.818 


3.540 
3.539 


6.82 
6.09 


Model Sn from Femandes & Monteiro (2003) 
this paper 


Model H- 




1.080 


1.667 


3.248 


5.79 


this paper 


Model H+ 




1.080 


2.007 


3.617 


6.33 


this paper 


Model A 


a Cen A 


1.111 


1.224 


1.5460 


6.96 


this paper 


M2000 A 
Th2002 A 


a Cen A 
a Cen A 
a CenB 


1.160 
1.100 
0.928 


1.228 
1.230 
0.867 


1.527 
1.519 
0.5025 


2.71 
4.85 
6.88 


Model Arv from Morel et al. (2000) 
Thevenin et al. (2002) 


Model B 


this paper 


M2000 B 
Th2002 B 


a CenB 
oCenB 
a CenB 


0.970 
0.907 
0.929 


0.909 
0.857 
0.869 


0.571 

0.5002 

0.4991 


2.71 
4.85 
5.86 


Model Brv from Morel et al. (2000) 
Thevenin et al. (2002) 


T2008 


Model M2 from Tans et al. (2008) 



TABLE 2 

Near-surface corrections 



Model 



(i^obs(n)) 
(mHz) 



{A!/obs(n)) 
(mHz) 



(t'i-cf(«)> 

(mHz) 



(mHz) 



a 
(mHz) 



Pre! 

(gem"'') 



(gem"'') 



Models 17-25 3034.15 
Models" 17-25 3034.15 
Model S"^ 17-25 3034.15 



Sun (1-0 = 3 100 /iHz) 

134.810 3038.95 135.854 
134.810 3386.40 150.761 
134.810 2540.67 114.155 



-4.73 1.00000 
-1.54 0.89644 
-8.67 1.19770 



1.408 
1.748 
0.982 



1.408 
1.405 
1.408 



H 


13-19 


1004.42 


57.244 


H 


13-19 


1004.42 


57.244 


)3 


13-19 


1004.42 


57.244 


H+ 


13-19 


1004.42 


57.244 



l3 Kyi (vo = 1000 pHz) 

1134.75 65.985 

1001.90 58.417 

948.58 55.714 

868.28 50.654 



-4.57 0.88978 

-5.19 1.00847 

-6.90 1.06723 

-5.32 1.16385 



0.329 
0.253 
0.226 
0.188 



0.260 
0.258 
0.258 
0.255 



Model S+ 16-26 2335.89 

M2000 A 16-26 2335.89 

Model A 16-26 2335.89 

Th2002 A 16-26 2335.89 



Q Cen A {uo = 2400 ^Hz) 

105.541 2509.96 114.026 

105.541 2362.88 107.949 

105.541 2345.37 106.344 

105.541 2303.36 104.927 



-3.18 0.93196 

-6.48 0.99141 

-2.05 0.99686 

-4.78 1.01626 



0.982 
0.882 
0.854 
0.833 



0.853 
0.867 
0.849 
0.860 



Th2002 B 19-27 3890.70 

Model B 19-32 4261.05 

T2008 19-32 4261.05 

M2000B 19-32 4261.05 

Model S" 19-32 4261.05 



a Cen B (i^o = 
161.482 3910.71 



161.988 
161.988 
161.988 
161.988 



4274.09 
4261.92 
4064.00 
3979.80 



4100 ^Hz) 

162.605 
162.913 
162.697 
155.400 
152.166 



-2.11 0.99534 

-1.87 0.99762 

-2.96 1.00087 

-4.15 1.05006 

-4.08 1.07225 



2.030 
2.009 
1.995 

1.819 
1.748 



2.011 
1.999 
1.998 
2.006 
2.010 
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Fig. 1 . — The difference between observed and calculated frequencies for radial modes in the Sun. The squares are for Model S, with the solid curve showing 
a fit to equation J4) with r = 1, which gives h = 4.90 (see Q. Also shown are the results of applying the same corrections to Model S" (crosses) and Model S+ 
(pluses). The dotted curves show the corrections calculated from equation (4). 
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Fig. 2. — The difference between observed and calculated frequencies for radial modes in /3 Hyi. The models shown are: Model H (squares). Model H 
(triangles). Model H+ (diamonds) and FM2003 (asterisks). The dotted curves show the corrections calculated from equation (4). 
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Fig. 3. — Same as Fig.|2] but for q Cen A. The models shown are: Model A (squares), M2000 A (triangles). Model S^ (pluses) and Th2002 A (asterisks). 
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Fig. 4. — Same as Fig.|2] but for a Cen B. The models shown are: Model B (squares), M2000 B (triangles). Model S" (crosses), Th2002 B (asterisks) and 
T2008 (diamonds). 



